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Abstract 

We consider two-dimensional percolation in the scaling limit close to criticality and use 
integrablc field theory to obtain universal predictions for the probability that at least one 
cluster crosses between opposite sides of a rectangle of sides much larger than the correlation 
length and for the mean number of such crossing clusters. 



1 Introduction 



Remarkable results have been obtained in the last two decades for crossing clusters in two- 
dimensional percolation at its critical point pc- In particular, following the numerical study of 
[U, Cardy used conformal field theory to derive an exact formula for the crossing probability 
Pj] that, in the continuum limit, at least one cluster spans between the horizontal sides of a 
rectangl^, a result later proved rigorously by other methods [5]. Cardy also determined the 
mean number of crossing clusters [6l [7]. 

For a rectangle of width L and height R, as a consequence of scale invariance both Py and 
Ny depend at pc only on the aspect ratio R/L. In this paper we consider these quantities in 
the scaling limit close to Pc, where, due to the presence of a finite correlation length ^ (much 
larger than the lattice spacing), they separately depend on L/^ and R/(,. We consider the limit 
L 3> 1^ and use boundary integrable field theory to determine the mean number of vertically 
crossing clusters, i.e. the clusters which span between the sides of the rectangle separated by 
the distance R, in the limit i? ^ ^. The result we obtain below pc is given in (|49|) . (|48|) . On 
the other hand, we can observe that for i? — )• oo below pc vertical crossing becomes extremely 
rare, so that P^, = Proh{Ny > 0) ~ Prob(A^^ = 1) ~ N^; from the leading term in (j48p we then 
obtain for the vertical crossing probability in the scaling limit below pc the universal result 

P,(L,i?)~A |e-«/«, R>L, (1) 

where 

A = ^i3-V3). (2) 

The correlation length ^ we refer to is defined by the decay of the probability P2{r) that two 
points separated by a distance r are in the same finite cluster: 

P2{r) r-" e-""/^ , r ^ oo , (3) 

with a = 1/2 below Pc and a = 2 above pc [8]; ^s related to the mass m appearing in (jl8|) and 
throughout the paper as 

J 1/m , p <Pc 



l/2m , p> Pc- 

We will also give a direct derivation of ([T]) which also yields the next term, given by (j5ip . in the 
large R expansion. The corresponding results in the scaling limit above pc are given in (j47p and 
(ESI). 

It has been previously known [Hi [TO] that for R = L ^ the crossing probability is a function 
of L\p—pc\'^ which decays exponentially to zero below pc and to one above, a feature investigated 
numerically in [TT| [T2| [TO] . 

We will start the analysis recalling in the next section the relation with the g-state Potts 
model and explaining how the latter is described in the scaling limit within the framework of 



^See [3] for the probability of simultaneous horizontal and vertical crossing. At pc results for the rectangular 
geometry are connected by conformal symmetry to other simply connected domains [2l [4] . 
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Figure 1: The different types of clusters whicli, depending on boundary conditions, determine 
the g-dependence of the Potts partition functions ([6117]) . 



boundary integrable field theory. In section 3 we take the limit g — )• 1 relevant for percolation 
and determine the quantities of our interest. 



2 Mapping to the Potts model and field theory 

It is well known that an efficient theoretical approach to random percolation is to see it as a 
limiting case of the g-state Potts model defined by the lattice Hamiltonian [14:\ [T5] 

n = -J^Ssi,sj , Si = l,...,q. (5) 

(id) 

The partition function X]{si} admits the expansion |16| — over bond 

configurations Q, where p = 1 — e~"^, n is the number of bonds in Q, n the complement to 
the total number of edges, and each cluster formed by adjacent bonds contributes a factor q 
corresponding to the number of color^ it can take; for g — )• 1 configurations are weighted as 
required for bond percolation. The above expression for the partition function holds as it is for 
free (f) boundary conditions. If instead the color of the spins on a boundary is fixed to be a, the 
clusters touching that boundary can take only the color a and do not contribute any factor q to 
the weight. In particular, if we denote by Z'^^ the Potts partition function on a rectangle with 
boundary conditions u, d, I and r on the upper, lower, left and right boundary, respectively, as 
already noted in [Bl [7] we have 

za = Y.p^'i^ - pf^'''^ = Y.p^^^ - (6) 

g g 
g g 



^Different values of the Potts spins can conveniently be associated to different colors. The Hamiltonian ([5]) is 
invariant under permutations of the colors. 
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where Nf, is the number of clusters which do not touch the horizontal boundaries, AT^ i^^d) the 
number of clusters which touch the upper (lower) but not the lower (upper) boundary, and Ny 
the number of those touching both horizontal boundaries (see Fig. [1]). It follows that the mean 
number of vertically crossing clusters can be written as 

7// 7// 

Ny = lim dq log 7f ff . (8) 

af fa 

Since boundary conditions on Potts spins loose physical meaning as g — t- 1 and sites do not 
interact in random percolation, ([8]) gives the mean number of clusters spanning between the 
horizontal sides of a rectangular window within the infinite plane on which the percolative 
transition actually takes place. 

Let us begin our field theoretical considerations for the scaling limit considering an infinitely 
long horizontal strip of height R. With imaginary time running upwards, the partition functions 
on the strip can be written as 

Zl^, = {B^\e-^^\B,)i^r, (9) 

where H is the Hamiltonian of the quantum system living in the infinite horizontal dimension, 
\Bd,u) are boundary states specifying initial and final conditions, and the vertical boundary 
conditions at infinity have the role of selecting the states which can propagate between the 
horizontal boundaries. Integrability of the scaling Potts model [T7] allows us to work in the 
framework of integrable field theories for which the bulk dynamics is entirely specified by the 
Faddeev-Zamolodchikov commutation rules (see e.g. |18| ) 



Alie,)A]ie2) = stfie, - 02)^],(^2)4(^i) > (lo) 

A\ei)A]i.02) = Sj'j{02 - 9i)A],{92)A'' (Oi) + 2^5(^1 - 62) , (11) 

where 4^(6') and A' (9) are creation and annihilation operators for a particle of species i with 
'f 



rapidit}o 0, and 5* ■-' {9) are two-body scattering amplitudes satisfying, in particular, unitarity 



Stfi0)sf:f{-9) = 6f5f (12) 

and crossing symmetry 

sff{9) = Sg(i7r-9). (13) 

Generic boundary states \Ba) can be written as superpositions of asymptotic states of the par- 
ticles created by A^, with vanishing total momentum in order to preserve horizontal translation 
invariance. The additional constraints coming from the requirement that a boundary condition 
preserves integrability were discovered in [19]. In particular, particles carrying momentum can 
only appear in pairs with vanishing total momentum. On the other hand, for the case of our 
interest of a theory satisfying 

Sff {0) = {-l) 6(6^ , (14) 



^Energy and momentum of a particle with mass m are given by (e,p) — (m cosh 6, msinh^ 
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states containing k > 2 particles of zero momentum are forbidden by (jlOp . Finally \Ba) takes 
the form 

\Ba)=Sae^p[- I —Va{e)\\n), (15) 



where \Q) is the vacuum state and 

5, = 1 + ^4(0), (16) 
Va{9) = K^i9)Ali-e)A]{0). (17) 

The boundary pair emission amplitudes Ka{0) satisfy equations involving the bulk amplitudes 
Sl^^ {9), and the constants follow from the relationj^ 

2iRese=oK'iO)=g:gi, (18) 

The exponential form of the boundary state is a consequence of the boundary Yang-Baxter 
equations, which give in particular [Va{d),Va{6')] = [Va{9),Sa] = [19] . 

We are now ready to use this formalism to evaluate the Potts partition functions entering 
([8]). The Potts field theory, i.e. the integrable field theory which describes the scaling limit of the 
Potts model in two dimensions, was solved exactly in [17] in the language of the spontaneously 
broken phase above Jc, in which the elementary excitations are kinks 4^(0)|Oq,) interpolating 
between degenerate ferromagnetic vacua \i^a) and |r2^) with different color. The vacua satisfy 
{Qal^is) = Sai3 and the admissible multi-kink states have the form 

4„+,a„(^n) . . . Ala,{02)Ai^a,m^a,) ■ (20) 

For these topological excitations the Faddeev-Zamolodchikov commutation relations take the 
form, 



AU0i)4'y(^2) = E - ^2)a1,(02)<(^i) > (21) 

s 

A^f,{0,)Al^{02) = Sj',{e2 - 9,)Alsi^2)As^{0i) + 2tt5^^5{9, - O^) , (22) 



where invariance of the theory under permutations of the colors allows for the four inequivalent 
scattering amplitudes represented in Fig. [21 they are given explicitly in [T7] and obey (1140 in 
the form 

5g(0) = (-1)5^5. (23) 

Both fixed (to a color a) and free (f) boundary conditions are integrable and the corre- 
sponding pair emission amplitudes were determined in [22]. They determine the boundary 
states \Ba) and \Bf) in the form that we now specify; in the following q will be parameterized 



*See I2ni[n] for the factor 1/2 in (0. 
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Figure 2: The four inequivalent kink-kink scattering amplitudes in Potts field theory (different 
indices denote different colors). 



as q = 4sin^(^), so that A — t- 1/2 corresponds to the percolation limit. For fixed boundary 
conditions we have 



with 



\Ba) = exp 



L2 



I a 



(24) 
(25) 



Kq{9) = itanh^-j exp 
nx{t) 



dt n\{t) 
t 2 cosh t 

sinh(| + ^)-sinh(f -^)_ 



sinh ( t 



2dt\ 
m J 



sinh(^) cosh(i) 



(26) 
(27) 



the integral in (|26p is convergent for 1/2 < A < 1. Kq{9) satisfies the boundary "cross-unitarity" 
relation 

Koi0) = K(20) + {q- 2)S^^Z{29)]Ko{-9) , (28) 

which together with (p3|) implies i^o(O) = 0, as already apparent from (p6|) : the absence of a 
pole at = explains 5^ = 1. 

For free boundary conditions we have instead 



1 



d9 



Vf(9) 



with 



vfi9) = y: 



K,{9)Alp{-9)A\^{9) + K2{9) 



(29) 



(30) 



5 



smh[A(4.^/3 - 26)] r(-4A/3 + 29 + l)r(7A/3 - 29) 

""'^'^ r(2A/3-2^ + l)r(A/3 + 2^) "^^'^ ' ^''^ 

sin^smh[A(OT-20)] sinh[A(4OT/3 - 20)] r(-4A/3 + 2^ + l)r(7A/3 - 2^) 
^~ sin^ smh(2A0) smh[X{-2iTT/3 + 29)] T{2X/3 - 29 + 1)T{X/3 + 29) 

(32) 

Q{9)=exp - I -- — — vii _2a; — VT 2Ay g-^^t^/^. ^ (33) 



where 9 = ^\ the integral in ([33p is again convergent for 1/2<A<1. In this case the residue 



at = is non-zero and give^ 

~g} = '-Rese=oK^{9) = ^Res,=oi^2(0) = (3 - g) ¥ r(l - |)r(|) ^^^^ ^^^^ 

2 2 4 A r(i + -]r(3) 

We also quote the boundary cross-unitarity conditions 

K,{9) = K(20) + {q- 2)S^^U29)]K,{-9), (36) 
K2{9) = [Sf^(20) + (g - 3)5f^(20)]E:2(-0). (37) 

3 Partition functions and final results 

In principle, the knowledge of the bulk and boundary amplitudes should allow the study of par- 
tition functions on the strip for any R through the boundary version [23] of the thermodynamic 
Bethe ansatz (TBA) [21]. In practice, however, the very non-trivial structure of Potts field 
theory seriously complicates the tasljfl. More pragmatically, here we plug the explicit expres- 
sions for \Ba) and \Bf) into ([9]) and exploit the fact that the states ([20|) are eigenstates of the 
Hamiltonian H with eigenvalues rn^^^^ cosh^j, m being the mass of the kinks. This leads to 
a large R expansion for the partition functions for which we compute below the terms coming 
from one- and two-kink states. 

Since we work in the kink basis, the partition functions we obtain in this way are those above 
Jc, that we denote Z^^'^, keeping the notation Z'^'^ for those below Jc- As an illustration, for z[a 
expansion of the boundary state leads to 

zLLiR) = 1 + I / -7:r^K{9')K,{9)e-^^''^°^'^' ^^^1^(9, 9') + ©(e"^-^) , (38) 

4 J-00 i^TTj ^^^^^ 

Map^^i9,9') = {n^\A^p{9')A^^{-9')Al^{-9)A\^{9)\n^) 

= (27r)2 [6(9' - 9)] + (27r)2 [8{9' + 9)f Sli{29') ■ (39) 



^There appears to be a typo in eq. (48) of |22j . In particular it does not reproduce §/(§) ~ 1 for the Ising 
model (q=2). 

®See |25l for the state of the art of TBA in the Potts model. 
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Figure 3: Diagrammatic interpretation of the two contributions entering the matrix element 



the last equality follows from formal use of (j22p and can be associated to the diagrams shown in 
Fig. [3l If L — )• oo denotes the horizontal size of the system, the squared delta functions in (j39p 
admit the usual regularizatiorlll 



[6(9' ±9)]- 



5{9' ±9)^ cosh 9; 



(40) 



free vertical boundary conditions have been imposed making no selection on the kink states 
which propagate between the horizontal boundaries. Exploiting boundary cross-unitarity 
and real analiticity Kq{—9) = Kq{9), G R, we then obtain 



Faa{R) 

Similarly one finds 



l + {q-l) mL Faa{R) + ©(e-^"^^) , 

oo 



mL > 1 , 







27r 



cosh 9 \Ko{9)\' e~ 



(41) 
(42) 



0(e 



^//(^) = - coshe (|iri(e)|2 + (g-2)|ir2(e)|2) e- 



2mR cosh ^ 



mL > 1 , (43) 
(44) 



zij{L,R) 
Faf{R) 



and Z 



ff 



Zj'-t. The partition functions Zf^^ with fixed vertical boundary conditions are 



l + {q-l)mLF^f{R) + 0{c 

— cosh9Re[K^{9)Ki{9)\ e" 
27r 



) , mL > 1 , 

2mR cosh 9 



(45) 
(46) 



obtained taking off from Z^^ the contribution of the states which are not of the form (|20p with 
ai = a and oin+i = /5- 



'It is important to stress that, as observed for other models in [23] (see also [26]), contributions to ((9} coming 
from states with more than two particles produce singularities whose regularization depends in general on the 
interaction. This is what makes difficult the determination of additional terms in the large R expansion within 
the approach we are following. 
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From these results we obtain for the mean number of crossing clusters in the scaling limit 
above Pc the universal result 

Ny{L, R) = lim dg log J/f ~ 1 + mL \<^{R) + 0(e-^"^)l , jnL > 1 , (47) 

^R) = Ae-™« + [Fff{R) + F^^{R) - 2F^f{R)]g=i , (48) 
where A = g'j\q=i reduces to ([2]). The additive term 1 in (I47p is produced by the overall factor 



q in (j43p and accounts for the contribution of the infinite cluster in the limit i? — )• oo. Notice 
that any normalization of the boundary states other than the one we used would anyway cancel 
in the combination of partition functions in ()47p . 

In order to determine the mean number of crossing clusters below pc we have to use the 
duality |14j of the Potts model to connect the partition functions (I6ll7|) below to the partition 
functions Z^^^ above Jc we have computed. Duality maps free boundary conditions into fixed 
boundary conditions and vice versa (see e.g. [27]). For our present purpose of counting the 
vertically crossing clusters, it is useful to observe that fixing the spins to the color a on both 
vertical sides, rather than leaving them free, has the only effect that the clusters touching at 
least one vertical side are not counted. Below pc, where all clusters are finite with a mean linear 
extension of order ^, such a boundary term does not affect Ny, which is extensive in L in the 
limit L/^ — 7- oo we are considering. So we can use ([8]) with the replacement // — t- aa in the 
vertical boundary conditions, and use dualitj|§ to obtain in the scaling limit below Pc 



NviL,R) 



yff yff 

lini dq log J' 
Zy Z^{ 

fa af J 



= mL [<^{R) + Oie-^""^)] , mL > 1 . (49) 

extensive part 

A different derivation of this result is given in the Appendix. 

The functions Fj f{R)\q=i and Faf{R)\q=i in (j^Hj) are well defined in spite of the poles at 
= in the amplitudes Ki[9) and K2{0) contained in the integrands in (I44p and ()46p . If the 
convergence of (I46p simply follows from Kq[Q) = 0, the case of Fjj is more subtle. Consider 
indeed the Laurent expansions Ki{e) = a-i/e + af + af e + i = 1, 2. If G M, due to 
the relation K*{9) = Ki{—9) the coefficients 02^-1 purely imaginary and the coefficients 
afl^ are instead real for all non-negative integers k; this in turn implies that the combination 
\Ki{0)\'^ — \K2{6)\^ entering Fff{R)\q=i does not contain any double or single pole at = 0. 

The function ()48p is plotted in Fig. HI where some numerical values are also listed. Since 
(j49p is extensive in L for any i?, it is tempting to check what our large R result gives in the 
conformal limit mR 0, for which the result iV^, ~ (V3/4)L/i? = (0.433. .)L/i?, L > i?, is 
known from [6]. Using the large 6 limits Kq — ;> e*'^/'^, Ki — ;> — 2e*''/^ and K2 — -y/Se*'"/^ at 
g = 1, (BSI) givefl 3/(27r)(L/i?) = (0.477.. )L/ii, with a 10% deviation from the exact result 
suggesting that (|48p may still provide a good approximation for mR of order 1. 



®In principle Z"" could be mapped into a linear combination of and a ^ ji. However, it follows from 
24| that the latter partition function vanishes identically on the infinitely long strip; at large L it is suppressed 
\s e 

^Given F{y) = dx coshx e^^""""^"^ f{x), with lim^j-^oo fix) = a, we have F{y) — >• a/y for y ^ 0. 
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Figure 4: Plot of the function ()48p ; few values are given in the table. 



We already explained how ([T]) follows from (I49p . The same result also follows from the 
observation that a lattice configuration with a vertical crossing is mapped onto a dual lattice 
configuration without horizontal crossings, and vice versa [15], so that|l^ 



P„ = 1 - Ph 



lim Z 



ff 



(50) 



The partition function Zjj is obtained picking up in (I43p only the contributions of the states 
compatible with the vertical boundary conditions a/3. In particular, since we are no longer 
summing over a and /3, the one-kink contribution in e~™^ now appears with multiplicity one 
rather than q(q — 1), and this gives ([I]) baclJ^. Concerning the two-kink contribution, only 
the term containing |i^2p survives now, again without the prefactor q{q — 1); |i^2p contains a 



singularity of the form 4g^ /9'^ at 6 = which has to be subtracted^i because produced by the 
propagation of states created by AQ^(0)Aj^^(0) which, as already observed, are not compatible 
with (^T]) and Hence, the contribution of order e~'^^l^ to be added to ([T]) is 



U{L,R) = -mL 



de 



cosh^ (|i^2(^)|g=ie 



-2m,Rcosh9 



4^2 



-2mR 



(51) 



The replacement 1/^^ — )• l/sinh^ 6, relevant for the convergence of the integral, comes from the 
fact that m J d9 cosh 9 = J dp, and in the momentum variable p the function \K2\'^ diverges as 

l/p2. 

The vertical crossing probability above pc is = 1 — lim^^^i Z^^. We already observed 
that Z^^ vanishes exponentially at large L, in agreement with the expectation that, due to the 
presence of an infinite cluster, above pc the crossing probability tends to 1 as we enlarge the 
window. More precisely, duality gives 



Pv{L,R) 



l-Ph{L,R) 
1 



A i^e-^/2« 



2^ 



U{R,L), 



L>R, 



(52) 
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In the continuum limit at pc (|50[) reproduces the known relation Pv + Ph = i [I] • 



Notice that our normalization of the boundary states ensures the conditions limq_>i 



required for percolation. 

^^The result obtained in [26] for simpler models (with purely transmissive scattering) by a TBA analysis amounts 
to such a subtraction. 
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Figure 5: The inset shows Monte Carlo data from |30j for the crossing probabihty below pc 
(diamonds) and for the complement to 1 of the crossing probability above pc (circles); the data 
refer to bond percolation on the square lattice with L = R = 256. The tails are plotted against 
mR using m = uiqIp — Pc\^^^, with niQ = 5.8. The continuous curve is the result ([1]), ([2]), (j52p . 
i.e. AmRe-""^. 

where the last line takes dU into account and holds at order e~^/^. 
Acknowledgments. We thank J. Cardy for a discussion about eq. (|5U|) . 

Note added. We learned from a referee that a scaling analysis of Monte Carlo data for Pv{L, R) 
in terms of a single scaling variable was performed in [30] and further discussed in |3H I32j. It 
is relevant for the present paper that the data of [30] allow a comparison with our results. The 
inset of Fig. [S] shows the data of [3D] for the crossing probability in bond percolation on the 
square lattice of size L = R = 256 lattice units; they satisfy the duality relation ([50|l for the 
crossing probability above and below pc (which for R = L specializes to P^, = 1 — P„) up to 
discrepancies to be ascribed to a mixing of finite size effects, corrections to scaling and statistical 
errors. In principle comparison of the data to ([T]), ([2]) and (|52p allows to fit the value of the only 
unknown parameter, i.e. the non-universal amplitude rriQ entering the relation m = ttiqIp — pd'^ , 
v = 4/3. We fit rriQ ~ 5.7 from the tail of the crossing probability below pc, and niQ ~ 5.9 from 
the tail above pc', consider that for R = 256 and niR around 10, ^ is around 25 below pc and 
around 12 above, so that the analysis is almost certainly affected by non-negligible corrections 
to scaling. On the other hand, the value 2,^2nd ~ ^■'^'^ measured in |33j for the amplitude of 
the second moment correlation length in bond percolation on the square lattice below pc, and 
it is known from [8] that mo = a,/i2nd ~ 5.4 a, with a equal 1 up to corrections that are not 
expected to exceed few percents. A comparison between the data for the tails and (HI, ([2]) is 
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given in Fig. \5\ the subleading term (j5ip is always very small and totally negligible in the range 
of mR shown in the figure. Putting all together, our conclusion is that the dat80 of [Ml are 
consistent with the results of this paper within the numerical uncertainties; an unambiguous 
verification will require simulations expressly targeting the tails on larger lattices. 

We are very grateful to H. Watanabe and C.-K. Hu for providing us with the data of j30j . 
and to the referee for bringing references [301 [3T| [32] to our attention and for noticing that the 

1 /2 

constant A, that we originally quoted in the form (|(2 — Vs)) , can equivalently be written 
as in ([2]). 



Appendix 

Let us denote by (j^abix) the field [2H] whose insertion at point x on the boundary changes the 
boundary condition from a to 6; of course (paa coincides with the identity /. If xi, . . . , 2:4 are 
the coordinates of the corners of the rectangle starting from the left upper corner and moving 
clockwise, we haveH 

Glurd = {4>luixi)(j)urix2)(prd{x3)(j)dliXi))j<j^ = Z''^^/ Z\\ , (53) 
Glurd = {4'lu{xi)cl)ur{x2)4>rd{x3)cl)dl{Xi))j*>J, = Z^udl " (^4) 

The fields (j)ab obey the natural operator product expansion [2] 

(t>af ■ (t>fl3 = Sal3 I + C{1 - 6al3)fJ'a(5 -\ , (55) 

that we write symbolically omitting the coordinate dependence, and using the notation fi^is = 
(j)ai3 for the kink field which switches between fixed boundary conditions with different colors. 
The field Ha/^ix) is dual to the Potts spin field o"a(x) = q6s(^x),a ~ 1) arid the relation 

{cra{x)crf^{y))j<j^ = {q5afs - 1) {n^s{x)fis-y{y)).J*>Jc (56) 

holds (see e.g. [29]). Since duality exchanges fixed and free boundary conditions, we then write 
the dual of as 

(f>fa- 4>af = I + C aa-\ , (57) 

with c' a new structure constant. For the boundary correlators Giurd, we have simple duality 
relations like Gfffa = Gaaaf, but also non-trivial ones like 

Gfafa = Ol Gafaf + ^2 Gafj3f , (58) 
Gfafl3 = as Gafaf + ^4 G^fpf , (59) 



Aspect ratios L/R > 1 were also analyzed in [30) . but for these cases the range of mL covered by the data is 
not large enough to allow comparison with (|52|l . 
^■^At Jc analogous relations were used in [2j[6]. 
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where a ^ j3. In order to determine the coefficients ai, . . . , 04 we use (j55p and (j57p to take the 
hmits xi — X2 and X3 — )• 3:4 on both sides, take into account to equate the coefficients of 
the two-point functions we are left with, and obtain 

01=03 = 1, a2 = {q-l){c /cf , a4^ = -{c/cf. (60) 

Since = 1 — hmg_^i zjj^ = 1 — hmq_j,i 03 + (^aZJj , comparison with ([50]) and ([60]) gives 
(d jc)^ = 1 at g = 1. Putting all together, the combination of partition functions in ([8]) can be 
written as 



'vffyff ^ t n 



Zff + 02 



(61) 



af fa ^J^'JJ^JJJ'^ ^ajaa^aaaj ^fa'^af 

with 02 = q — 1 + 0{{q — l)"^) . Now it is not difficult to use our expressions for the partition 
'^""^ to check that limq_>.i dq ] 



functions Z\\ to check that lim„_>.i d„\ogR gives the r.h.s. of ([19]) . 
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